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Using lattice quantum chromodynamics (QCD) we perform an unquenched calculation of the 
gluon propagator in Landau gauge. We use configurations generated with the AsqTad quark action 
by the MILC collaboration for the dynamical quarks and compare the gluon propagator of quenched 
QCD {i.e., the pure Yang-Mills gluon propagator) with that of 2+1 flavor QCD. The effects of the 
dynamical quarks are clearly visible and lead to a significant reduction of the nonperturbative 
infrared enhancement relative to the quenched case. 
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I. INTRODUCTION 



The gluon propagator, the most basic quantity of 
QCD, has been subject to much calculation and spec- 
ulation since the origin of the theory. In particular there 
has long been interest in the infrared behavior of the Lan- 
dau gauge gluon propagator as a probe into the mecha- 
nism of confinement 111. Some authors have argued it to 
be infrared finite 0| while others favored infrared 
singular U [(J. There is a longhistory of its study on 
the lattice, in quenched QCD 0, H |1, El E1MJI1 
Q E E [H [H EH and in quenched SU(2) [U Hj. 
The restriction to quenched lattice gauge theory calcu- 
lations has been due to the lack of sufficient computa- 
tional resources. The quenched theory differs from full 
QCD only in the relative weighting of the background 
gauge configurations (due to the fermion determinant), 
but the evaluation of the Green's functions is otherwise 
the same. In the quenched approximation the fermion 
determinant is replaced by unity and this corresponds 
to the complete suppression of all quark loops. The re- 
moval of quark loops is equivalent to the limit where all 
sea-quark masses are taken to infinity. In this paper, we 
report the first results for the gluon propagator from an 
unquenched lattice computation. 



II. DETAILS OF THE CALCULATION 



The gluon propagator is gauge dependent and we work 
in the Landau gauge for ease of comparison with other 
studies. It is also the simplest covariant gauge to imple- 
ment on the lattice. Landau gauge is a smooth gauge 
that preserves the Lorentz invariance of the theory, so it 
is a popular choice. It will be interesting to repeat this 
calculation for the Gribov-copy free Laplacian gauge, but 
that will be left for a future study. 

The MILC configurations were generated with the 
0(a 2 ) one-loop Symanzik improved [24[ Liischer-Weisz 
gauge action |25| . The dynamical configurations use the 
"AsqTad" quark action, an 0{a 2 ) Symanzik improved 
staggered fermion action. (3 and the bare sea-quark 
masses are matched such that the lattice spacing is held 
constant. The lattices we consider all have the same di- 
mensions. This means that all systematics are fixed; the 
only variable is the addition of quark loops. The param- 
eters are summarized in Table [I] The lattice spacing is 
approximately 0.125 fm [26| . 

In Landau gauge the gluon propagator is entirely trans- 
verse. In Euclidean space, in the continuum, the gluon 



We study the gluon propagator in Landau gauge u sing 
configurations generated by the MILC collaboration |22| 
available from the Gauge Connection (2?|]. These use 
"AsqTad" improved staggered quarks, giving us access 
to relatively light sea quarks. We find that the addition 
of dynamical quarks preserves the qualitative features of 
the gluon dressing function q 2 D(q 2 ) in the quenched case 
- enhancement for intermediate infrared momenta fol- 
lowed by suppression in the deep infrared - but produces 
a clearly visible effect. A significant suppression of the 
infrared enhancement with respect to the quenched case 
is observed. It is interesting to compare these results to 
those of a recent Dyson- Schwinger equation study p3| . 



TABLE I: Lattice parameters used in this study. The dy- 
namical configurations each have two degenerate light quarks 
(up/down) and a heavier quark (strange). In physical units 
the bare masses range from ~ 16 MeV to ~ 79 MeV. The 
lattice spacing is a ~ 0.125 fm. 

Dimensions (5 Bare Quark Mass # Configurations 
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x 64 


8.00 


quenched 


192 
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20 :i 


x 64 


6.76 


0.01, 0.05 


193 
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20 3 


x 64 


6.79 


0.02, 0.05 


249 
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20 3 


x 64 


6.81 


0.03, 0.05 


212 
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20 :i 


x 64 


6.83 


0.04, 0.05 


337 



2 



propagator has the tensor structure 



D nv(q) = ( <V q2 



and at tree-level, 



(i) 



(2) 



With this lattice gauge action the propagator at tree-level 
is 
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where 
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a is the lattice spacing and is the length of the lattice 
in the /x direction. As explained in Ref . Q > this suggests 
a "kinematic" choice of momentum, 
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(5) 



ensuring that the lattice gluon propagator has the correct 
tree-level behavior. 

The bare gluon propagator, D{q) is related to the 
renormalized propagator Dn(q;n) through 



D(q) = Z 3 ((j,,a) D R {q;n) 



(6) 



where fi is the renormalization point. In a renormalizable 
theory such as QCD, renormalized quantities become in- 
dependent of the regularization parameter in the limit 
where it is removed. Z% is then defined by some renor- 
malization prescription. We choose the momentum space 
subtraction (MOM) scheme where Z3 (//, a) is determined 
by imposing the renormalization condition 



1 

7 2 



(7) 



i.e. it takes the tree-level value at the renormalization 
point. In the following figures we have chosen = 4 
GeV. 



III. SIMULATIONS RESULTS 

Lattice studies strongly suggest that the quenched 
gluon propagator is infrared finite [l4|. As is custom- 
ary, we will begin by considering the (necessarily finite) 
gluon dressing function, q 2 D(q 2 ). In Fig. ^ we compare 
the well-known quenched dressing function with that for 
2+1 flavor QCD. For the moment we only consider the 
lightest of our dynamical quarks as we expect that they 
will show the greatest difference from the quenched case. 
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FIG. 1: Gluon dressing function in Landau gauge. Full trian- 
gles correspond to the quenched calculation, while open circles 
correspond to 2+1 flavor QCD. As the lattice spacing and vol- 
ume are the same, the difference between the two results is 
entirely due to the presence of quark loops. The renormaliza- 
tion point is at /i = 4 GeV. Data has been cylinder cut |lql . 
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FIG. 2: The sea-quark mass dependence of the Landau gauge 
gluon propagator dressing function renormalized at = 4 
GeV. Filled squares correspond to u and d bare masses ~ 63 
MeV and bare s-quark mass ~ 79 MeV. Open circles corre- 
spond to the same strange-quark mass, but with bare u and d 
masses ~ 16 MeV. Data has been cylinder cut |l6f . Increasing 
the sea-quark masses alters the results in the expected way, 
i.e. towards the quenched data. 



Indeed there is a clear difference between quenched and 
dynamical quark behavior in the infrared region. The ad- 
dition of quark loops to the gluon propagator softens the 
infrared enhancement without altering its basic features. 
The screening of dynamical sea quarks brings the 2 + 1 
flavor results significantly closer to the tree-level form, 
q 2 D{q 2 ) = 1. 

In Fig. we show the gluon dressing function for the 
lightest and for the heaviest u and d quark masses in our 
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FIG. 3: The renormalized propagator at one momentum point 
in the infrared hump of the gluon dressing function (q ~ 1.12 
GeV) is shown here as a function of the bare light-quark mass. 
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FIG. 4: The sea-quark mass dependence of the Landau gauge 
gluon propagator renormalized at /i = 4 GeV. Filled trian- 
gles illustrate the quenched propagator while filled squares 
correspond to bare up/down masses ~ 63 and bare strange- 
quark mass ~ 79 MeV. Open circles correspond to lighter 
bare up/down masses ~ 16 MeV but with the same strange 
quark mass. Data has been cylinder cut [T^. 
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FIG. 5: The light sea-quark mass dependence of the renormal- 
ized gluon propagator at a momentum point in the infrared 
region (q ~ 0.31 GeV). 

point from the infrared hump (q ~ 1.12 GeV) and plot 
it for each choice of bare light-quark mass. Although the 
variation in the propagator at this momentum is only 
4.5% over the range of quark masses investigated here, 
the light sea-quark mass dependence is clearly resolved. 

In Fig. 0] we present results for the gluon propagator, 
D(q 2 ). The largest effects of unquenching are observed in 
the deep infrared. The shape of the curves suggest that 
previous results indicating the infrared-finite nature of 
the quenched gluon propagator [J^ are unchanged upon 
unquenching. The results suggest that the gluon propa- 
gator of QCD is infrared finite. It will be interesting to 
examine the behavior of -D(O) as a function of volume to 
elucidate this aspect of the gluon propagator further. 

Finally, in Fig. [5] the light sea-quark mass dependence 
of the renormalized gluon propagator is illustrated for a 
momentum point in the infrared region. To avoid finite 
volume artifacts, the second smallest nontrivial momen- 
tum is considered. Whereas the mass dependence of the 
propagator for the masses studied here is at the 4.5% level 
for q ~ 1.12 GeV, the variance is larger in the infrared 
region at 6% for q ~ 0.31 GeV. 



set. These correspond to bare light-quark masses of ~ 16 
MeV and ~ 63 MeV respectively; a factor of four differ- 
ence. The bare strange-quark mass is the same in both 
cases (~ 79 MeV). The mass dependence of the gluon 
dressing function is only just detectable. We expect that 
increasing the sea-quark masses further will interpolate 
between the curves in Fig. ^ We see that the gluon prop- 
agator changes in the expected way. As the sea-quark 
mass increases, the curve moves toward the quenched re- 
sult. However, for the range of bare quark masses studied 
here the change is relatively small. This transition would 
be better studied with heavier sea quarks. 

Another view of the mass dependence of the gluon 
propagator is provided in Fig. [3] We choose one data 



IV. CONCLUSIONS 

The addition of quark loops has a clear, quantitative 
effect on the gluon propagator. While its basic structure 
is qualitatively similar there is significant screening of the 
propagator in the infrared. As anticipated, the effect is 
to suppress the non-abelian enhancement of the gluon 
propagator in the nonperturbative infrared-momentum 
region. This is relevant to analytic studies of the gluon 
propagator and confinement 23]. Despite the clear dif- 
ference between the quenched and dynamical results, we 
see little dependence on the dynamical quark mass for 
the range of available light sea-quark masses. The depen- 
dence that is observed is consistent with expectations. 
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Calculations on finer lattices arc currently being made, 
which will provide more information on the ultraviolet 
nature of the propagator and provide a test for finite 
lattice spacing artifacts. We would like to extend the 
study to a wider range of dynamical masses to study 
both the chiral limit and the transition to the quenched 
limit. Finally, a study of the volume dependence of the 
propagator will provide valuable insights into the nature 
of the propagator at q 2 = 0. 
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